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Abstract 

We show that the density /U of the Smith normal form (SNF) of a random integer matrix exists and equals 
a product of densities //p® of SNF over Z/p^Z with p a prime and s some positive integer. Our approach 
is to connect the SNF of a matrix with the greatest common divisors (gcds) of certain polynomials of 
matrix entries, and develop the theory of multi-gcd distribution of polynomial values at a random integer 
vector. We also derive a formula for pps and compute the density p for several interesting types of sets. 
Finally, we determine the maximum and minimum of pps and establish its monotonicity properties and 
limiting behaviors. 


1. Introduction 

Let M be a nonzero nxm matrix over a commutative ring R (with identity), and r be the rank 
of M. If there exist invertible n x n and m x m matrices P and Q such that the product PMQ 
is a diagonal matrix with diagonal entries di, d 2 , ■ ■ ■ ,dr,0,0,... ,0 satisfying that di \ dj+i for all 
1 < i < r — 1, then PMQ is the Smith normal form (SNF) of M. In general, the SNF does not 
exist. It does exist when i? is a principal ideal ring, i.e., a ring (not necessarily an integral domain) 
for which every ideal is principal. This class of rings includes the integers Z and their quotients 
Z/gZ, which are the rings of interest to us here. In fact, for the rings Z/gZ we will be particularly 
concerned with the case g = p®, a prime power. For principal ideal rings, the diagonal entries are 
uniquely determined (up to multiplication by a unit) by Pi-idi = Qi {1 < i < r), where go = 1 and Pi 
is the greatest common divisor (gcd) of aWixi minors of M. We have the following correspondence 
between the SNF and the cokernel of M: coker M ~ R/diR © R/d 2 R © • • • © R/d^R © R^~^. 

There has been a huge amount of research on eigenvalues of random matrices over a held (see, 
eg., n. 0. |10j . [I2]). Less attention has been paid to the SNF of a random matrix over a 
principal ideal ring (or more general rings for which SNF always exists). Some basic results in this 
area are known, but they appear in papers not focused on SNF per se. We develop the theory 
in a systematic way, collecting previous work in this area, sometimes with simplihed proofs, and 
providing some new results. 

We shall dehne the density /i of SNF of a random nxm integer matrix as the limit (if exists) 
as fc —)■ CX3 of the density of SNF of a random nxm matrix with entries independent and 
uniformly distributed over {—k,—k + l,...,k} (see Dehnition 13.11 below for a precise dehnition). 

As a motivating example, the probability that di = 1 for a random nxm integer matrix is 
the probability that the nm matrix entries are relatively prime, or equivalently, that nm random 
integers are relatively prime, and thus equals l/({nm), where C(') is the Riemann zeta function. 

If we regard the minors of an n x m matrix as polynomials of the nm matrix entries with integer 
coefficients, then the SNF of a matrix is uniquely determined by the gcds of the values of these 
polynomials (recall the dehnition of SNF from the beginning). This inspires us to study the theory 
of multi-gcd distribution of polynomial values. 
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Given a collection of relatively prime polynomials in Z[a:i, X 2 , ■ ■ ■, Xd], let g{x) be the gcd of the 
values of these polynomials at a; = {xi,X 2 , ■ ■ ■ ,Xd)- We shall dehne the density A of g{x) of a 
random d-dimensional integer vector x as the limit (if exists) as fc —)■ cxd of the density of g{x) 
with X uniformly distributed over {—k, —k -f 1,..., k}'^ (see Dehnition l2.1l for a precise dehnition). 

In the spirit of previous work in number theory such as 0 . la. ca and the Cohen-Lenstra 
heuristics (la. 0 ). one might conjecture that A exists and equals the product of density Ap of g{x) 
over {h/phY over all primes p. In fact, we will prove this conjecture with the more general density 
Xps of g{x) over 'L/p^'L for sets of form (I2.5p (see Theorem 12.81) . with the aid of a result in number 
theory m Lemma 21]. Note that the special case that s = 0 or 1 follows from P, Theorem 2.3] 
directly. In particular, this result applies to the probability that g{x) = 1, in other words, that 
the polynomial values are relatively prime. Furthermore, all these results hold for the multi-gcd 
distribution of polynomial values, namely, when g{x) is a vector whose components are the gcds 
of the values of given collections of polynomials at x. 

Then we apply this theory to the SNF distribution of a random integer matrix to show that the 
density p (of SNF of a random n x m integer matrix) equals a product of some densities Pps of 
SNF over Z/p^Z for sets of form fl3.4|) (Theorem 13.81) . We also derive a formula for pps (Theorem 
13.2p . which allows us to compute Pps and hence p explicitly (Theorem 14.3|) . Some special cases of 
this formula coincide with [161 Exercise 1.192(b)] and P pp.233, 236]. Another paper related to 
our work is [HI- 

On the strength of these results, we determine the value of p for some interesting types of 
sets, specihcally, matrices with hrst few diagonal entries given, matrices with diagonal entries all 
equal to 1, and square matrices with at most i{= 1, 2,..., n) diagonal entries not equal to 1, i.e., 
whose corresponding cokernel has at most £ generators; further, for the last set we establish the 
asymptotics of p as £ —)■ 00 . In the case of £ = 1 (which is equivalent to the matrix having a cyclic 
cokernel), our results echo those of Ekedahl [61 Section 3] via a different approach. We also show 
that the probability that a random integer matrix is full rank is 1, and that p of a hnite set is 0. 

Additionally, we hnd the maximum and minimum of pps{D) over all diagonal matrices D ] 
whereas regarding it as a function of p,s,m,n and D, we hnd its monotonicity properties and 
limiting behaviors. 

The remainder of this paper is organized as follows. Section [2] develops the theory of multi- 
gcd distribution of polynomial values. Section [3] applies this theory to the SNF distribution and 
derives a formula for pps. Section 0] computes the density p for several types of sets. Finally, 
Section 0] determines the maximum and minimum of Pps and discusses its monotonicity properties 
and limiting behaviors. 

We shall assume that throughout this paper, p represents a prime, pj is the j-th smallest prime, 
and Hp means a product over all primes p. 

2. Multi-gcd Distribution of Polynomial Values 

Suppose that d and h are positive integers and Fi, F 2 ,..., G Z[a;i, X 2 ,..., x^] are nonzero 
polynomials. Let 

g{x) := gcd(Fi(a;), F 2 (t), • • ■ , Fh{x )), x G Z'^ 

be the gcd of the values Ei(x), F 2 {x),, Fh{x), and g{x) = 0 if Fj{x) = 0 for all I < j < h. 

We shall dehne the density of g{x) of a random d-dimensional integer vector x as the limit (if 
exists) of the density of g{x) with x uniformly distributed over {—k, —k -|- 1,..., := Z^^^ as 

/c —)■ cxD, precisely as follows. 

Definition 2.1. (i) For Z C Z, we denote by X^^\Z) the probability that g{x) G Z with x 
uniformly distributed over Z^^^. If lim/j^oo X^’^^Z) = X{Z) exists, then we say that the probability 
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that g{x) G Z with x a random d-dimensional integer vector is X{Z). If this is the case, then 
X{Z) E [0,1] since X^^l{Z) E [0,1] for all k. 

(ii) We dehne similarly the gcd distribution over the ring of integers mod p^: for prime p 
and positive integer s, we denote by X^pJ{Z) the probability that g{x) E Z (mod p^) (up to 
multiplication by a unit) with x uniformly distributed over and by Xps(Z) the probability 
that g(x) E Z (modp'^) (up to multiplication by a unit) with x uniformly distributed over ('L/p^'LY. 

More generally, for a hnite set V of prime and positive integer pairs (p, s) (with p a prime and 
s a positive integer), we denote 

Pv-= n p' 

{p,s)€'P 

and by A^(Z) the probability that g{x) E Z (mod Pp) (up to multiplication by a unit) with x 
uniformly distributed over Z^^^, and by Xpp(Z) the probability that g(x) E Z (mod Pp) (up to 
multiplication by a unit) with x uniformly distributed over (Z/PtcZ)'^. Note that Xpp(Z) is the 
number of solutions to g(x) E Z (mod Pp) (up to multiplication by a unit) divided by The 
situation discussed in the previous paragraph is the special case that V consists of only one element 
(p, s) and Pp = p^. 

(hi) The above dehnitions also extend to the distribution of multi-gcds. Suppose that U = 
{Ui}^=i is a collection of w nonempty subsets Ui of {Pi, P 2 ,..., Fh}- Let 

(2.1) gi(x) := gcd (F(x) : F E Ui), x G Z'^ 
and 

g(x) := (pi,P 2 , • • ■,9w)(x) E IP, 

then we adopt the above dehnitions of functions X^^^ , A, Ap^ and Ap^ for Z with only one 

slight modihcation: replace “up to multiplication by a unit” with “up to multiplication of the 
components of g by units”. 

For convenience, we shall always assume that the notion g(x) G Z (mod Pp) implies the eguiv- 
alence of multiplication of its components by units and that the random vector x is uniformly 
distributed on its range (if known, e.g., or (Z/PpZ)'^). 

Remark 2.2. The density Ap(-) dehned above in Dehnition l2.1l fiii is consistent with the normalized 
Haar measure on Z^ , as in [T^ . 

In this section, we establish the properties of Ap^ and A, the existence of A, and a connection 
between A and the Aps’s. Then we apply these results to determine the probability that the 
polynomial values are relatively prime. 

2.1. Multi-gcd Distribution over Z/PpZ. 

We show that the density A^(-) over Z^^^ (dehned above in Dehnition 12.11) converges to the 
density App(-) over Z/PpZ as fc —)■ 00 , and that App(-) equals n(ps)e'P 

Theorem 2.3. For any Z C Z*^, we have 

(2.2) Ap,(2)= 5 ; Ap,({4) 

z^Z (mod P-p) 

lim aW(Z) = Ap,(Z)= n V(U. 

(p,s)G'P 


and 

(2.3) 
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Proof. (1) The first equality fl2.2p follows directly from Definition 12.11 

(2) For the second equality of fl2.3p . we let Np.p{Z) be the number of x G (Z/DpZ)^ for which 
g{x) G Z (mod Pp). The Chinese remainder theorem along with Dehnition 12.11 then gives 

P^\p^{Z) = Np^{Z)= J] Nps{Z)= J] = H Xps{Z). 

{p,s)gP {p,s)gP {p,s)&P 

Dividing both sides by Pp leads to the desired equality. 

(3) For the hrst equality of (j2.3|), we hrst observe that if p | 2/c + 1, then X^p\z) = Ap,(^) 
by dehnition. If p f 2fc + 1, then we proceed by approximating 2/c + 1 by a multiple of Pp and 
estimating A^(Z) using App(Z). 

Let A; G Z such that K ■.= 2k + 1> Pp, then there exists q G Z+ such that 
(2.4) q ■ Pp < K < {q + 1) ■ Pp . 


It follows that for any integer p, there are either g or g + 1 numbers among Z(fc) that equal y mod 
Pp. Thus the number of x G Z(),^ for which for which g(x) G Z (mod Pp) is between g'^iV' and 
(g + where N' := Np.p{Z), therefore 

-qdN' (g + l)'^iV'' 




J^d 


Rd 


. Jh 


Thanks to fl2.4p . we have 

Jk c 


- 1 

{q + l)^N'' 


\( ^ 

fq + lYN'' 

[{q + l)Pvt 

{qPv)^' _ 


\q + l) P$'^ 

1 _ 


whose left and right endpoints both converge to N'/Pp as g —>■ cxo. Hence 

X%(Z) ^ N' !P$ = Xp„(Z), 
as g ^ CX3 , or equivalently, as /c —)■ cx3 , as desired. 


□ 


2.2. Multi-gcd Distribution over Z . 

We show some properties of the density A of set unions, subtractions and complements. They 
will be very useful in determining the value of A for specihc sets (such as in Remark 12. 9l fiiiH. 

Theorem 2.4. Suppose that {Za}a£A (d^re pairwise disjoint subsets of IP such that A(Z„) exists 
for all a E A. If A is a finite set, then 

^ (Cag^ Za) = ^ ^ X(ZA) ■ 
a^A 

Proof. By Dehnition 12.11 we have 

VA(Z„) = V hm X^^^Zp,,) = hm = hm A^"^) (U«eA 

' ^ ^ k^oo k^oo * ^ k^oo 

q:Ev4 aGA Oi^A 

and the conclusion follows. □ 

Theorem 2.5. Suppose that Z' O Z P IP such that X{Z') and X{Z) both exist, then 

X{Z \ Z') = X{Z) - X{Z'). 

In particular, for the complement Z^ of Z in IP, we have 

A(Z") = 1 - X{Z). 
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Proof. By Definition 12.11 we have 

X{Z) - X{Z') = lim - lini X^’^^Z') = lim = lim X^^\Z \ Z') 


k—yoo 

and the conclusion follows. 


/c—>-oo 


k—^oo 


□ 


Theorem 2.6. Suppose that y such that A(3^) = 0, then for any Z <Py, we have X{Z) = 0 
as well. 

Proof. Since X^’^^Z) > 0, Z C y and linifc^.oo = A(3^) = 0, we obtain 

0 < liminf A(^)(Z) < limsup A(^)(Z) < limsup A(^)(y) = A(y) = 0. 

fc—>oo fc—>oo 


Therefore 
as desired. 


lim X^’‘\Z) = 0, 

k^oo 


□ 


2.3. Connection between A and Xps. 

We show that the density A exists and in fact, equals the product of some Apds. 

Assumption 2.7. For all 1 < i < tc, we have 

gcd(Fi, F 2 ,..., Fh) = gcd (F : F e Ui) = 1 in Q[xi, X 2 , ...,Xd]. 

Theorem 2.8. Suppose that Assumption \2. 7| holds. Given positive integers r <w and yi, I <i < 
r, let y = Y^=i'P^j Pj ^he j-th smallest prime and Sj nonnegative integers, j = 1,2,... such 
that yi I y for all 1 < i < r, then the probability X{Z) exists for 

(2.5) Z = {{zi, Z 2 ,..., Zn,) e'Ll : Zi = yi, \/i < r] , 
and in fact 

00 

(2.6) A(Z) = n\->«(Z). 

i=i ' 

Remark 2.9. (i) The right-hand side of fl2.6l) is well-dehned since Aps(-) e [0,1] for all p and s. 

(ii) The special case that all s/s are either 0 or 1 follows from [6l Theorem 2.3]. 

(hi) We have assumed that the i/i’s are positive. In fact, in the case that i/j = 0 for some i, we 
have A(2^) = 0 on the strength of Theorem 12.61 and that the probability that a nonzero polynomial 
at a random integer vector equals zero is 0 (see Theorem 12. 15l lii) below). 

To prove Theorem 12.81 we need Theorem 12.31 and the following two lemmas. 

Lemma 2.10. ([TU Lemma 5.1] or [151 Lemma 21]) Suppose that F,G E L[xi, X 2 , ■ ■ ■, Xd] are 
relatively prime as elements of Q[xi,X 2 ,... ,Xd]. Let be the probability that p \ F{x), G{x) for 
some prime p > i with x uniformly distributed over i.e., 

:= # {x e Lfi^^ : 3 prime p > i s.t. p \ F(x), G(x)} / {2k + 1)'^. 

Then 

(k) 

lim lim sup = 0 . 

fc—AOO 

Lemma 2.11. Suppose that Gi,... ,Gh E Q[xi,X 2 ,...,x^] {h> 2) are relatively prime, then there 
exists V = {vs,..., Vh) E L^~‘^ such that 

h 

Gi,G 2 + ViGi 

i=3 



gcd 


1 . 
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Proof. We prove by induction on h. The case h = 2 is trivial since gcd(Gi, G 2 ) = 1- 
Base case: h = 3 . 

We prove by contradiction. Assume the contrary that 

gcd(Gi, G 2 + z Gf) 7 ^ 1, V z G Z . 

Suppose that the polynomial factorization of Gi is 0 i 02 ''' 4>u, then each G 2 + z G^ is a. multiple 
of some factor (j)u{z) of Gi (1 < u{z) < u). Since there are inhnitely many 2 ;’s, by the pigeonhole 
principle, at least two of the u^zfs are the same, say u{z) = u{z') {z 7 ^ z'). Then 

4>u{z) I {G 2 + 2 G 3 ) — (G 2 + z'Gs) = (z — z') G 3 

thus (fu(z) I G 3 and hence (/>u(z) I (^*2 + 2 ; G 3 ) — ZG 3 = G 2 • Recall that 4>u(z) \ Gi as well. This 
contradicts with the condition that Gi,G 2 and G 3 are relatively prime. 

Inductive step-, from h — 1 to h (> 4). Assume that the statement holds for h — 1. 

Let H -.= (G 2 , G 3 ,..., Gh) and Hi := Gi/H (2 <i<h), then 

(2.7) gcd(Gi, H) = gcd(Gi, G 2 , • • • , G;,) = 1 = gcd(i 72 , • • • , . 

According to the induction hypothesis for H 2 , H 3 , • • • , if/j, there exists v = (^ 4 ,..., Vh) G 
such that 

h 

( 2 . 8 ) H',-.= H3 + J 2 ^iH, 

i=A 

satishes 

gcd(i72,i^3) = l- 

Combining with fl2.7p gives 

gcd{G 3 ,G 2 ,H',H) = gcd{G 3 ,H 2 H,H',H) = gcd{G 3 ,gcd{H 2 H,H’,H)) = gcd{G 3 ,H) = 1 . 

Thus we can apply the base case h = 3 to Gi, G 2 , H'^H to get an integer 2 ; such that 

gcd(Gi,G 2 +^77^77) = !. 

Finally, we represent H'^H back to a linear combination of the Gfs with integer coefficients by 
dehnition fl 2 . 8 p : 

h h 

77'77 = 77377 + Y,v^H,H = G 3 + J^P.G,, 

2=4 2=4 

therefore 

gcd ^Gi, G 2 + ^ G 3 + z PiGj^ = 1 , 

namely, the statement holds for h with the new v = {z, ZV 4 ,..., zvh). □ 

Now we are ready to prove Theorem 12.81 
Proof of Theorem \2.8\. Let 

•= {(Pi) + l)}i=i 5 7 e Z+ , 

then Theorem 12.31 gives 

t 

K(2)=nw‘G)' 

i=i ' 
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Since A G [0,1] for all j, we can let £ ^ cx): 


(2.9) 


lim Ap,, (2) = n K-‘*' (2) = RHS ot JZg ■ 

)-00 ^ J-J- p,. 


Therefore it snffices to show that 

( 2 . 10 ) 


i=i 


lim A<*l(2) = lim \p (Z ). 

/e—>-oo >-oo ^ 


Since i/ is hnite, there exists j* G such that Sj = 0 for all J > f- Let 

X Z2j • • • y ZiD ) G Z_j_ .2(1 ^2 ■ ■ ■ 1} 

then for any j > j*, the sets Z and X are equivalent mod pj under multiplication of the components 
by units. 

We dehne for £ > j*, 

A{t) := {x G : g{x) G Z (mod XpJ} , := {x G : g{x) G Z (mod XpJ} , 

:= {x € Zf,y. g(x) e Z} (c T‘)(«)) , 

and 

:= AlX)(£) 

then 

( 2 . 11 ) 


X^^\Z) = 


Rd 


W(u = 




with K ■= 2k + 1. Therefore 


( 2 . 12 ) 


(^) = lim Xp! (Z) = lim 


Rd J^d 


^ fc—AOO '' ' fc—AOO 

Combining with the hrst equation in fl2.1ip leads to 


Rd 


lim sup 

fc—>-oo 


A.q3),lims„pMh±|X!T).Ap IZ) 


k^oo 


and 


Rd 

Uminf A<‘-(^) > limmf - Um.up^ > Ap, (^) -lim.up 

fc^oo An-oo A“ i. .np, ^ i. A“ 


Once we show that 
(2.13) 


hmhmsnp^^O, 


taking £ —)■ cx) in the above two inequalities will yield fl 2 . 10 p . 

Now let us prove (I2.13p . For any x G , there exists j > £ {> j*) such that g{x) ^ X (mod 

= Pj) (recall that Z and X are equivalent). Hence pj \ gr^{x) for some p < r. 

Recall that grf is the gcd of some relatively prime Fj’s. If two or more Fj’s are involved, then 
applying Lemma 12.111 to these Fj’s leads to two relatively prime linear combinations and of 
these FiA with integer coefficients. If there is only one Fi involved, then it must be a constant since 
the gcd of itself is 1 in Q[xi, X 2 ,..., x^]. In this case, we take = Fj so that gcd(^^, Rri) = 1 

still holds. 

Since pj \ gn{x), we have pj | Qri{x),'Hn{x). Hence 

r 

w, 


B“>(f)C U 


(2.14) 
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where 


:= [x e : 3 j > ^ s.t. Pj \ ^^(x) ,. 
Applying Lemma 12.101 to and Hr] gives 

(£) 

lim lim sup-— = 0 , V n . 

fc^oo K<i 


Combining with fl2.14p . we obtain 


lim sup lim sup 

£^oo k^oo K 


< lim sup lim sup - 

i^oo k^oo 1 7 £—>-oo fc^oo 

77=1 77=1 


y7 1 ™ 1 ™ sup 


#<V) 

Rd 


= 0 . 


On the other hand, since > 0, we have 


hminfhmsup^^^!^>0. 

•^—>■00 J^^QQ l\ 


Hence fl2.13p indeed holds. 


□ 


2.4. Relatively Prime Polynomial Values. 

An interesting application of Theorem 12.81 is to determine the probability that the polynomial 
values are relatively prime. 


Theorem 2.12. Let w = 1 and Ui = {{Fi,F 2 , ■ • • iFh}} in Definition \2. li 

(a) //Pi, F 2 ,... ,Fh are not relatively prime in Q[xi, X 2 ,..., x^], then A({1}) = 0 ; 

(b) If Fi, F 2 ,..., Fh are relatively prime, i.e., 

(2.15) gcd(Pi,P 2 , • • • ,Ph) = 1 in Q[xi,X 2 , ... ,Xd]. 


then we have 

(i) A({1}) exists and 

A({i})=nMpHi 

p 

(ii) the asymptotic result 

(2.16) A,({0}) = 0(p-^); 


(hi) A({1}) = 0 if and only if Xp{{l}) = 0 for some prime p, i.e., if and only if there exists a prime 
p such that Pi(x), F 2 {x),..., Fh{x) are multiples of p for all x ; in words, the probability that the 
values of relatively prime polynomials at a random integer are relatively prime is 0 if and only if 
there exists a prime p such that these polynomials are all always multiples of p. 


Remark 2.13. Theorem I2.12l fbhiii and Lemma 12.141 in the proof below are special cases of the 
Lang-Weil bound mi Theorem 1]. We present a considerably simpler and more approachable 
proof. As mentioned in Remark 12.91 and [T^ Remark of Lemma 21], Theorem l2.12l fbRil follows 
from [ 6 l Theorem 2.3]; whereas its special case h = 2 was shown in [HI Theorem 3.1]. 

Proof, (a) Let G = gcd(Pi, P 2 , ■ ■ ■, Ph), then G is a non-constant polynomial. If the gcd g{x) = 1, 
then G(x) = ±1. Thus A*^^)({1}) < + <^g=-i 0 as / —)■ CX 3 on the strength of Theorem 

I2.15l fii). where (c = ±1) is the probability that G(x) = c with x uniformly distributed over 
Hence A({1}) = 0 . 

(b) (i) follows directly from Theorem 12.81 For (ii), we prove by induction on d. First, we notice 
the following facts: 

1. If h = 1, then Pi must be a constant due to Assumption fl2.15l) . Hence Ap({0}) = 0 for all 
p > I Pi I and fl2.16l) follows. 
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2. If /i > 2, by Lemma 12.111 there exist two linear combinations G and H of the Fj’s with 
integer coefficients snch that gcd{G,H) = 1 in Q[a:i,X 2 ,...,xj. Then p\g{x) implies that 
p I gcd(G(x), H{x)), so it snffices to prove for the case h = 2. 

3. Assnme that h = 2. Let L be the greatest total degree of the F^s. If L = 0, then Fi, F 2 and 
thns g are nonzero constants. Thus <Jp = 0 for any p > \g\ and fl2.16p follows, so we only need to 
prove for L > 1 . 

Base case: d = 1. Assume that h = 2 . 

Thanks to Assumption fl2.15p . there exist Hi,H 2 G Z[xi] such that HiFi + F12F2 = G with G 
a positive integer constant. If p \ g{x), then p | C as well. Hence = 0 for all p > C and (12.161) 
follows. 

Inductive step: from d — 1 to d (> 2 ). Assume that the statement holds for d — 1 and that h = 2 
and L > 1 . 

Since L > 1 , without loss of generality, we can assume that Fi is not a constant and xi appears 
in Fi. We recast Fi as a univariate polynomial Gi G (Z[x 2 ,..., Xrf])[xi] of degree Li, i = then 
Li > 1. Let G Z[x 2 , • • •, (0 < j < Li) be the coefficients of x\ in Gi. 

Since Fi and F 2 are relatively prime in Q[xi, X 2 ,..., x^] by Assumption (12.151) . we have 

gcd ( 7 ij :l<f< 2 , 0 <j<Lj) = l = gcd(Gi, G 2 ) in (Q[x 2 ,..., Xd])[xi] . 

As a result, there exist Hi, H 2 G (Z[x 2 ,..., X(i])[xi] and ids ^ Z[x 2 ,..., xj such that HiGi + 
H 2 G 2 = ids and (idi, H 2 , ids) = 1 in Q[xi, X 2 , ...,Xd\. 

If p I p(x), then p I {Gi{x 2 , ..., x^)) (xi) (V i), ids(x 2 ,..., x^) and either 

( 1 ) p I 7 jj(x 2 ,..., Xd) for all i and j; or 

(2) p f 7 ij(x 2 , ...,Xd) for some i,j. 

Case (1). Recall that Li > 1. By the induction hypothesis for the at least two polynomials: 
7 i,j(x 2 , • • •, Xd) (1 < i < 2,0 < j < Li), the probability that Case (1) happens with (x 2 ,..., Xd) 
uniformly distributed on (Z/pZ)*^”^ is 0 (p“^). 

Case ( 2 ). We need the following asymptotic result. 

Lemma 2.14. Let G G Z[xi, X 2 , ■ ■ ■, Xd\ he a nonzero polynomial, p a prime, and ap the probability 
that p I G(x) with x uniformly distributed over ifLlplL)'^, then we have 

(2.17) ap = 0{p-^). 

Proof. Let L be the total degree of G. If L = 0, then G is a nonzero constant. For any prime 
p > G, we have Up = 0, thus (I2.17P holds. 

Now we assume that L > 1. We prove by induction on d. 

Base case: d = 1. 

Since the number of roots of G mod p is at most L, we get ap < L/p and hence (I2.17p . 

Inductive step: from d — 1 to d (> 2 ). Assume that the statement holds for d — 1. 

We recast G as a univariate polynomial Gi G (Z[x 2 , X 3 ,..., Xd])[xi]. Let 71 G Z[x 2 ..., xj be the 
leading coefficient of Gi. Observe that the total degree of Gi is at most L. If 71 (x 2 ,..., Xd) ^ 0 
(mod p), then the probability that p|Gi(xi) with xi uniformly distributed over Z/pZ is no 
greater than L/p, according to the base case d = 1. On the other hand, the probability that 
P I 71 ( 2 ^ 2 , • • • 5 Xd) with (x 2 ,..., Xd) uniformly distributed over (Z/pZ)'^“^ is 0 (p“^) by the induction 
hypothesis for 71 . Combining these two cases, we conclude that the probability that p | G(xi, X 2 , ■ ■ ■, 
Xd) with (xi,X 2 ,... ,Xd) uniformly distributed over (Z/pZ)'^ is at most L/p + 0(p“^) = 0(p“^). 
Hence the statement holds for d, as desired. □ 
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Now we go back to the proof of Theorem l2.12l fbViib Thanks to Lemma 12.141 the probabil¬ 
ity that p I ■ ■ ■ ,Xd) with (x 2 ,... ,Xd) uniformly distributed on is more¬ 

over, for each {x 2 , ■ ■ ■ ,Xd) that satishes p f 7 ij(T 2 ,...,x^) for some i,j, the probability that 
p I (Gi(x 2 ,..., Xrf)) (xi) with xi uniformly distributed on Z/pZ is 0{p~^). Hence the probabil¬ 
ity that Case (2) happens with (xi,X 2 , ... ,Xd) uniformly distributed on (Z/pZ)'^ is (0(p“^))^ = 

0 {p-^). 

Combining Cases (1) and ( 2 ), we conclude that the statement holds for d as well, as desired, 
(hi) If Ap({l}) = 0 for some prime p, then A({1}) = 0 by (i). 

Now assume that Ap({l}) > 0 for all prime p. On the strength of (ii), there exist a positive 
constant c and a positive integer j* such that 


Thus 


Pj* > 1 + 0 y/c, Xpj ({0}) <CPj'^ , V j > j*. 


n (1 - Y({o})) > n 




3=3 


1-4 

p] 


> 1 


E C 

-2 > 1 

P) 


EA 


3=3 


i=Pi 




> 0 , 


where in the second inequality, we take advantage of the well-known inequality: 


(2.18) (1 - ^i)(l - ^2) • • • (1 - > 1 - - ^2- 

for 61 , 62 ,... , 6 u e [0)1]) which can be proved easily by induction on u (base cases: u = 1,2; 
inductive step from u to n -|- 1: (1 — 5i)(l — ^ 2 ) ■ ■ ■ (1 — (5„+i) > (1 — (5i)(l — (52 — ■ • ■ — <5m+i) > 
1 — (5i — (^2 — • • • — (5u+i). 

Hence 

j*-i 00 

=n Y({i}) ■ n (1 - y({o})) > o- 


2 . 5 . Zero Polynomial Values. 

Remark I2.9l fiiii used a well-known result that the probability that a nonzero polynomial at a 
random integer vector equals zero is 0 ([HI Lemma 4.1]). We conclude this section with a different 
proof by estimating this probability from above by ap and applying Lemma 12.141 

Theorem 2.15. Let G G Z[xi,X 2 ,... ,Xd\ be a nonzero polynomial, p a prime, ap^'^ the probability 
that p I G(x) with x uniformly distributed over Z^^^, and ap the probability that p | G(x) with x 
uniformly distributed over {fL/p'LY, then 

(i) we have 

a^Y'^ ap as k ^ 00 , and a^^'^ < 2^ap , V /c > (p — l)/2 ; 

(ii) the probability a^^'^ that G{x) = 0 with x uniformly distributed over Z^^^ goes to 0 as k ^ 00 ; 
in words, the probability that a nonzero polynomial at a random integer vector eguals zero is 0. As 
a conseguence, for any given integer c, the probability that G{x) = c is either 0 or 1 (consider the 
polynomial G{x) — c). 

Proof, (i) We follow a similar approach as in the proof of the hrst equality of fl2.3p . Let /c G Z 
such that K ■.= 2k + 1 > p. Then there exists q G Z+ such that 

(2.19) q-Pv <K <{q + l)-P-p. 
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<J. 


{k) 


Kd 


( 2 . 20 ) 

Thanks to (I2.19p . we have 

(g + lYp’^ap 


Rd 


:= S, 


( 2 . 21 ) 


Sfc C 


_[iq + l)pY' {qpY 


q + l 


a. 


p > 


g + 1 


(Tri 


whose left and right endpoints both converge to CTp as g —)■ cxd . Hence 


—)■ cr. 


p ’ 


as g ^ oo , or eqnivalently, as fc —)■ oo 


It follows that for any integer y, there are either g or g + 1 nnmbers among that eqnal y mod 
p. Farther, the nnmber of x G (Z/pZ)'^ for which p\G{x) is p'^Cp , thus the number of x G Z(),^ for 
which p I G{x) is between q^p'^cTp and (g + l)'^p'^cTp . Therefore 

q’^p’^ap (g + lYp^a.^ 


Additionally, we deduce < 2'^crp from fl2.20p and fl2.2ip along with g > 1. 

(ii) The probability is no greater than <jY\ which by virtue of (i) and Lemma [2. 141 converges 
to 0 as p, A; —)• oo with k > {p — l)/2 . □ 


3. SNF Distribution 

Let m < n he two positive integers. We shall define the density of SNF of a random n x m 
integer matrix as the limit (if exists) of the density of SNF of a random nx m matrix with entries 
independent and uniformly distributed over Z(fc) as A: —)■ oo (see Dehnition 13.11 below for a precise 
dehnition). 

If we regard the minors of an n x m matrix as polynomials of the nm matrix entries with integer 
coefficients, then the SNF of a matrix is uniquely determined by the values of these polynomials. 
Specifically, let Xi, X2 ,..., Xnm be the nm entries of an n x m matrix, F)’s be the minors of an n x m 
matrix as elements in Z[xi, X2 ,..., Xnm], Ui be the set of i x i minors (1 < i < m), then the SNF 
of this matrix is the diagonal matrix whose A-th diagonal entry is 0 if gi{x) = 0 and gi(x)/gj_i(x) 
otherwise, where x = (xi, X2 ,..., Xnm) and gi{x) is defined in fl2.ip . 

In this spirit, the multi-gcd distribution as well as the results in Sections l2.1112. 3l have analogues 
for the SNF distribution of a random integer matrix. This section presents these analogues and 
the next section will use them to compute the density p for some interesting types of sets. 

Conventionally, the SNF is only defined for a nonzero matrix; however, for convenience, we shall 
define the SNF of a zero matrix to be itself, so that SNF is well-defined for all matrices. This 
definition does not change the density (if exists) of SNF of a random n x m integer matrix since 
the probability of a zero matrix with entries from Z(fc) is 1/(2A; -|- which converges to 0 as 

k ^ oo. 

We denote the SNF of an n x m matrix M by SNF(M) = (SNF(M)jj)n,xm and let S be the 
set of all candidates for SNF of an n x m integer matrix, i.e., the set of n x m diagonal matrices 
whose diagonal entries (di, d 2 , ■ ■ ■, dm) are nonnegative integers such that dj+i is a multiple of d*, 
i = 1, 2,..., m — 1. 

For ease of notation, we shall always assume that the matrix entries are independent and uni¬ 
formly distributed on its range (if known, e.g., Z(fc) or 'LjP'p'L)-, and that the notion SNF(M) G iS or 
SNF(M) = D (mod PY) for some 5 C §, D G § and P^p = n(ps)e'P^'* ^ the equivalence 

of multiplication of the entries of M by units in Z/P-pZ, thus we can assume for convenience that 
the entries of SNF(M) (mod Pp) are zero or divisors of Pp. 

Definition 3.1. (i) For S C §, we denote by the probability that SNF(M) G S with 

entries of M from Z(fc). If lim^^oo/i(^^(iS) = p{S) exists, then we say that the probability that 
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SNF(M) G S with M a random nxm integer matrix is /i(iS). If this is the case, then /i(iS) G [0,1] 
since G [0,1] for all k. 

(ii) We dehne similarly the SNF distribution over the ring of integers mod p^: for prime p and 
positive integer s, we denote by p^pJ{S) the probability that the SNF(M) G S (mod p*) with 
entries of M from and by Pps{S) the probability that SNF(M) G S (mod p®) with entries of 
M from Z/p*Z. 

More generally, for a hnite set V of prime and positive integer pairs (p, s) (with p a prime and s 
a positive integer), we denote by fip^{S) the probability that SNF(M) G S (mod Pp) with entries 
of M from Z(fc), and by fipp{S) the probability that SNF(M) G S (mod P-p) with entries of M from 
Z/P-pZ. Note that ppp{S) is the number of matrices M over Pp such that SNF(M) G S (mod Pp) 
divided by Pp’". The situation discussed in the previous paragraph is the special case that V 
consists of only one element (p, s) and Pp = p^. 


In this section, we establish a formula for /Ups, discuss the properties of /up^ and /i, show the 
existence of p and represent it as a product of Ppds. 

3.1. SNF Distribution over Z/PpZ. 

We have the following formula for jUps and analogue of Theorem 12.31 for SNFs. 


Theorem 3.2. (i) Given a prime p, a positive integer s and a seguence of integers 0 = oq < ui < 
■■■< Os < Os+i = m, let a := (ai, 02 ,..., Us) and Da E S be the diagonal matrix with exactly 
(at — Oj-i) p*“^ ’s, i.e., a* non-p^-multiples, 1 < i < s on its diagonal. Then we have 


(3.1) 

where 


Tp-{{Da]) =P 




[p,n][p, m] 


[p, n - Os] [p, m - a^] HLi \Py “ ®*-i] ’ 


[p,o] = i, [p,^] := JJ (1-p ^) , iez+. 

j=i 

(ii) For any S PE), we have 


TPv (‘^) 

DgS (mod P-p) 


and 

(3.2) 


lim = 1‘PriS) = n 1‘p’iS ). 

{p,s)&P 


Proof, (ii) and (hi) are direct applications of Theorem 12.31 to SNFs. For (i), we compute the 
number of n x m matrices over Z/p^Z whose SNF is Da by [3 Theorem 1] (or |S1 Theorem 2]) 
and simplify it to 


(3.3) 

Thus 


pELi [{n+m)ai 



n!:(;n;=r“‘(i-p-d 


=: N. 


TA{Da]) = p-^"^N = RHS of ([ST]) . 


□ 


Remark 3.3. In the case of s = 1, the formula fl3.3p gives the number of n x m matrices over 
Z/pZ of rank oi and is consistent with [161 Exercise 1.192(b)]; whereas in the case of u = m, a 
calculation shows that fl3.3p is consistent with the results in P pp. 233, 236] (their lAutif] is our 
N). 
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3.2. SNF Distribution over Z . 

The properties of A of set unions, subtractions and complements in Section 12.21 also carry over 
to SNFs. They will be useful in determining the value of fi for some specihc sets (for instance, the 
singleton set of the identity matrix as in Section 03]) ■ 

Theorem 3.4. Suppose that are pairwise disjoint subsets of § such that fi{Sa) exists for 

all a & A . If A is a finite set, then 

aeA 

Theorem 3.5. Suppose that 5' C 5 C § such that and /i(iS) both exist, then 

^,{S \ S') = /i(5) - /x(5'). 

In particular for the complement S'^ of S in §, we have 

f.iS") = 1 - /i(5). 

Theorem 3.6. Suppose that T G S such that fi(T) = 0, then for any S C T, we also have 

M5) = o. 

3.3. Connection between pi and pips. 

Theorem 12.81 has an analogue for SNFs as well, by virtue of the following well-known lemma (see 
[21 Theorem 61.1] for an easy proof). 

Lemma 3.7. Fix a positive integer r. The determinant of an r x r matrix as a polynomial of its 
entries Xi,X 2 , ■ ■ ■, Xr 2 is irreducible in Q[a;i, X2 ,..., x^i] ■ 

For any i < m A {n — 1) (i.e., min {m,n — 1}, which is m if m < n, and n — 1 if m = n, 
recalling that m < n), the set f/j contains at least two different minors, which are both irreducible 
as polynomials of the entries on the strength of Lemma 13.71 and therefore relatively prime. Hence 
Assumption 12.71 holds with w = m A {n — 1). This allows us to apply Theorem 12.81 to SNFs and 
obtain the following analogue. In addition, we will compute the density p^{S) explicitly later in 
Section 14.11 

Theorem 3.8. Given positive integers r < m A {n — 1) and di 1^2 | • • • | dr, let z = YY^=i'P^j 
Pj the j-th smallest prime and Sj nonnegative integers, j = 1,2,... such that dr\z, then the 
probability pi{S) exists for 

(3.4) S = {D := (Aj)nxm e S : A,i = d*, V i < r} , 
and in fact 

OO 

(3.5) /x(5) = JJ/i .,.+i(5). 

.=1 ' 

Remark 3.9. (i) The right-hand side of (13.51) is well-dehned since p.ps{-) G [0,1] for all p and s. 

(ii) We have assumed that r < m A (u — 1); in fact, we have /x(iS) = 0 otherwise. Recall that 
m < n and note that r < m, thus in the case of r > m A (u — 1), we must have r = m = n. 
As a result, any matrix M with SNF(M) G S satishes \M\ = ±A. We will show later that the 
probability that the determinant of a random n x n integer matrix equals ±c is 0 for all constant 
c (Theorem 14.5p . 

(hi) We have also assumed that the dfs are positive; in fact, we have pi{S) = 0 otherwise. If 
dj = 0 for some i, then all i x i minors of any matrix M with SNF(M) G S are zero. Applying 
Theorem 14.51 to c = 0 yields the desired result. 
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4. Applications 

Now we apply Theorems 13.21 and 13.81 to compute the density /i explicitly for the following subsets 
of §: matrices with hrst few diagonal entries given (i.e., with the form of fl3.4p h full rank matrices, 
a hnite subset, matrices with diagonal entries all equal to 1, and square matrices with at most 
i{= 1,2,... ,n) diagonal entries not equal to 1. 

4.1. Density of the Set 03.41) . 

For the set S of 03.4p . i.e., of matrices with hrst r diagonal entries given, we take z = dr in 
Theorem 13.81 then it suffices to compute /ips+i(iS) for each (p, s) = {pj,Sj). In mod the set 
S has m — r + 1 elements (see 04.121) below). Further, since formula 03.1|) gives the density Pps+i 
of each element of S, one can take the sum over S to get an expression for /ips+i(iS) (Theorem 
13.2p . and compute this sum explicitly when m — r is small, such as in Theorems 14.81 and 14.91 below. 
However, this sum is hard to compute when m — r is large, for example, when m is large and r is 
hxed; in this case, we recast S as the difference between a subset of § and the union of other r — 1 
subsets such that for each of these r sets, its density Pps+i is given directly by 03.ip . 

We work out two examples to illustrate this idea, and then deal with the general case. 


4.1.1. The First Example: Relatively Prime Entries. 

Our approach reproduces the following result mentioned at the beginning of this paper. 


Theorem 4.1. Let S he the set of 03.4p with r = 1 and di = 1, then we have 


(4.1) 


r{S) 


1 

({nm) ’ 


where (('(•) is the Riemann zeta function. 


Proof. Applying Theorem 13.81 with r = 1, di = 1 and Sj = 0, j = 1,2,... gives 

(4.2) = 

P 


therefore it reduces to computing /ip(iS) for each p. 

Recall the equivalence of multiplication by units, therefore we only have two choices for matrix 
entries in mod p: 1 and 0. The set S (mod p) consists of all the matrices in § whose hrst 
diagonal entry is 1, thus S = {Da '■ o, = (oi) > 1} (mod p) (recall from Theorem 13.21 that 
a := (oi, 02 ,..., ttg) and that Da € S is the diagonal matrix with exactly Oj non-p®-multiples on 
its diagonal). Therefore 

/4p(‘5) 1 /ip({Z)(g)}) . 

We apply fl3.1|) to get pp({D(o)}) = thus fip{S) = 1 — p”"""*. Plugging into fl4.2p along with 

the Euler product formula 


(4.3) 


n(l-P'‘)= Ae(O.l). Vi>2 

P ^ V / 


yields fHTl) . 


□ 


4.1.2. Another Example. 

Theorem 4.2. Let S he the set of 03.41) with r = 2, di =2 and d 2 = 6, then we have 

(4.4) /i(5) = P22 (5) /i32 (*5) JJ /ip(*S), 

p>3 
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where 

(4.5) /X22(5) = 2-""* (^1 - 2-^^ - ^2-~^ 

(4.6) fi 32 {S) = (1 _ 3-(n-i)(m-i)^ ^)(l^-3 

1 3 


(4.7) fip{S) = 1 - - p-(—1)(—1) . P P 


1 - 


(n-l)m nm-l 

Y, P-+ Y p~‘' 


(n—l)(m—1) n(m—1)+1 


Proof. The first equation fl4.4p follows directly from Theorem 13.81 with r = 2,di = 2 ,d 2 = z = Q, 
Si = S 2 = 1, Sj = 0, j > 3. Therefore it reduces to calculating Ppo{S) for {p,s) = (2,2), (3,2) 
and (p, 1) with p > 3 . 

Case 1. p > 3 and s = 1. 

Recall the equivalence of multiplication by units, therefore we only have two choices for matrix 
entries in mod p\ 1 and 0. The set S (mod p) consists of all the matrices in S whose hrst two 
diagonal entries are 1, thus S = {Da '■ cl = (cq) > 2} (mod p) (recall a again from Theorem I3.2p . 
Therefore 

hp(‘5) = 1 — /rp({Tl(o)}) — /rp({Tl(i)}). 

We then apply (13.1 p to get /ip({i7(o)}) and Pp({i7(i)}), and (14.7p follows. 

Case 2 . p = 2 and s = 2 . 

We have three choices for matrix entries in mod 2^: 1, 2 and 0. The set S (mod 2^) consists of 
all the matrices in S whose hrst two diagonal entries are 2 , thus S = {Da={ai,a 2 ) : Oi = 0 , 02 > 2 } 
(mod 2^). Therefore 

(4-8) ^- 22 ( 5 ) = /I22({71(ai,a2) : Ol = 0}) — /r 22 ({71(0,0)}) — h 22 ({-D( 0 ,l)}) • 

Notice that the set {D(^ai,a 2 ) : cq = 0} (mod 2^) consists of all the matrices in S whose diagonal 
entries are all multiples of 2 (i.e., either 2 or 0 ); in other words, in mod 2 , it contains only one 
element - the zero matrix. Hence 


h22({-D(ai,a2) : — 0}) — P 2 {{D(o)}) . 

Plugging into (|4.8|) and applying (|3.ip to get /i2({71(o)}), h22({71(o,o)}) and p22({71(o,i)}), we obtain 

(O. 

Case 3. p = 3 and s = 2. 

We have three choices for matrix entries in mod 3^: 1, 3 and 0. The set S (mod 3^) consists of all 
the matrices in § whose hrst two diagonal entries are 1 and 3, respectively, thus S = {Da={ai,a 2 ) ■ 
oi = 1, 02 > 2} (mod 3^). Therefore 

(4-9) h32(‘5) = P32{{D(^ai,a2) ■ = 1}) “ h32 ({-D(l,l)}) . 

Notice that the set {D(^ai,a 2 ) : cti = 1} (mod 3^) consists of all the matrices in S whose hrst diagonal 
entry is 1 and all other diagonal entries are multiples of 3 (i.e., either 3 or 0); in other words, in 
mod 3, it contains only one element -- the diagonal matrix whose diagonal entries are 1, 0, 0,..., 0. 
Hence 

P'32({77(ai,a2) • 1}) h'3({7)(l)}) • 

Plugging into (14.9p and applying (13.ip to get /i3({i7(i)}) and /i 32 ({Zi)(i i)}), we obtain (14.61) . □ 
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4.1.3. The General Case. 


Theorem 4.3. Let S be the set of fl3.4p in Theorem \3.8\ with dr = YlJLiP’jC then for {p,s) = 
{pj, Sj), j = 1,2,..., we have 


Pps+-L{S) = p 




\p,n][p, m] 


\p,n- a^] [p, m - a^] HLi \P^ “ ®*-i] 


r—1 


(4.10) 


p-G-i)im-e)-Y.Ui{n-ai){m-ai) _ 


\p, n] \p, m] 


[p,n- i][p, m- i][p,i- Os] ni=i[P, Oj - Oi-i] ’ 


where Oj (0 < i < s) is the number of non-p'^-multiples among di, d 2 , ■ ■ ■ ,dr (thus Os < r — 1) . In 
particular, when s = 0 (which holds for all but finitely many j ’s), we have 


(4.11) 


/^j(5) = 1 - 

1=0 


\p,n — T\ [p, m — i] [p, £] 


The value of p{S) is then given by Theorem \3.8\ with z = dr . 


Proof. Recalling from Theorem 13.21 the notation of Da, we recast S as 

(4.12) S = {Da=(ai,a 2 ,...,as+i) : o* = Oj , 1 < f < s , Os+i > r} (mod , 
and therefore 

r—1 

(4.13) fipB+ll^S) Pps+l ({T^a.=(ai,a 2 ,...,as+l) ' ’ 1 ^ ^ s}) ^ ^ ({-^(ai,Q 2 i....Qs,^)}) • 

l=as 

Notice that the set {Da={ai,a 2 ,...,as+i) • > 1 < * < -s} (modp®"''^) in the hrst term on the right- 

hand side of fl4.13p consists of all the matrices in § with exactly a, (1 < i < s) non-p*-multiples on 
its diagonal; in other words, in mod p^, it contains only one element ~ the diagonal matrix with 
exactly di non-p*-multiples, i.e., (di — di-i) p*“^’s, 1 < z < s on its diagonal. Hence 

(4.14) Mp‘’+^({^Da=(ai,a 2 ,...,as+i) • 0 ,i di , 1 ^ z < s}) hp® ({- 0 (ai,a 2 ,...,as)}) • 

Plugging into (| 4.13 p and a pplying ([3T]) to get /ip<>({T>(ai,a 2 ,...,a,)}) and /ip»+i({T>(ai,a 2 ,1 < 
£ < r — 1, we obtain fl4.10p . □ 

Remark 4.4. We notice that the density Pps({Da}) of fl3.ll) is a polynomial of p~^ with integer 
coefficients since m — Og + X]i=i(®« ~ Oi-i) = The Pp(S) of fl4.11l) is also a polynomial of p~^ 
with integer coefficients and with constant term 1 (see the Pp{S) of 04.71) as an example). If we 
replace each occurrence of p by where z is a complex variable, and plug into 03.5p of Theorem 
13.81 we get an Euler product for some kind of generalized zeta function. 

For instance, when m = n = 3, for the set S in Theorem 14.21 we apply 04.7|) to get 

PpiS) = 1 - p“^ - p~^ - p“® -h p“^ + p“® = (l - P~^) (l - P~^) (l + + P~^) ■ 

Taking the product over all primes p and applying the Euler product formula 04.3p . we obtain 

Plugging into 04.4p . we see that to obtain the density p(S), it reduces to computing np(l + 
p~^), or to understanding the Euler product np(l + P~^^ -\-p~^^). 

It would be interesting to study whether such an Euler product for some generalized zeta function 
(1) has any interesting properties relevant to SNF; (2) extends to a meromorphic function on all 
of C ; (3) satisfies a functional equation. 

































SMITH NORMAL FORM DISTRIBUTION OF A RANDOM INTEGER MATRIX 


17 


4.2. The Determinant. 

The determinant of an m x m matrix can be regarded as a polynomial G of its entries. Note 
that G is not a constant since it takes values 1 and 0 for the identity matrix and the zero matrix, 
respectively. Thus we can apply Theorem 12.151 to G and obtain the following. 

Theorem 4.5. Let c be an integer. The probability that the determinant equals c for an m x m 
matrix with entries from Z(fc) goes to 0 as k ^ oo; in other words, the density of the determinant 
of a random mx m integer matrix is always 0 . 

This result plays an important role in the next two theorems. The hrst of them shows that the 
probability that a random n x m integer matrix is full rank is 1 . 

Theorem 4.6. If S S satisfies Dm,m = 0 for all D = {Dij)nxm £ S, then we have /i(iS) = 0 ; 
in other words, the probability that an n x m matrix with entries from Z(fc) is full rank goes to 1 
as k ^ oo . 


Proof. If SNF(M)m,m = 0 , then all m x m minors of M are zero. Therefore the result follows from 
Theorem 14.51 with c = 0 . □ 


When m = n, we can generalize Theorem 14.61 to S with hnitely many values of Dm,m’s. 


Theorem 4.7. Suppose that m = n and iS C §, then we have fi{S) = 0 if the set : D = 

{Dij)nxn £ *5} is finite; in particular, this holds for any finite subset 5 C § . 

Proof. For any M such that SNF(M) = D E S, we have \M\ = ±Di iD2^2 ■ ■ 'Dn,n- As a conse¬ 
quence, if Dn^n = 0, then \M\ = 0; if Dn^n 7^ 0, then the Di/s are divisors of and therefore 
\M\ has hnitely many choices. The result then follows from Theorem 14.51 □ 


If Dn^n 7 ^ 0 for all D E S, then we have another proof of Theorem 14.71 without invoking Theorem 
14.51 We cannot take advantage of fl3.2l) from Theorem 13.81 since r = m = n>mA{n — l)m this 
case. Instead, we will start from the observation that with P{i) a product of 

primes and / the identity matrix, then bound from above by 2 "^/ip(£)({/}) through a 

similar idea as in the proof of fl2.3p (approximating 2fc -|- 1 by a multiple of P{i)), and hnally show 
that /ip(^)({/}) —)■ 0 as £ —)■ cxo. 


Another Proof of Theorem Jf..! with 7 ^ 0 for all D E S . Let / be the n x n identity matrix 
and j* E Z_|_ such that pj > c for all j > j*. Then for any j > j*, SNF(M) G S (mod pj) only if 
SNF(M) = / {mod Pj). 

Applying fl3.3l) with s = 1 and ai = n (or [161 Exercise 1.192(b)]), we get the number of n x n 
non-singular matrices over 'Ljpj'L'. 

pf\pj,n] := fij. 

Set 


P{i) ■.= Pj*Pj*+i---pe, I>j*. 

Then SNF(M) G 5 (mod P{i)) only if SNF(M) = / (mod P{i)). Hence < lx^P^f^{{I}). 

By the Chinese remainder theorem, the number of n x u matrices over 'LlP{P)'L whose SNF 
equals / mod P{i) is 


(4.15) 




I I 

j=j* j=j* 


For any integer k with K ■.= 2k + 1 > P{i), there exists q E such that 
(4.16) q-P{i) < K < (q + l) ■ P{i). 
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Then for any integer z, there are at most q + 1 numbers among that equal z mod P{t)- 
Therefore the number of n x n matrices with entries from Z(fc) whose SNF is equal to I mod P{i) 
is at most (g + Hence 


(4.17) 




(q + ifm 


{q+ir"m 

[qP{e)f 


q + l 


m 

p{iy- 


< 2" 




n\ 


j=r 


on the strength of fl4.16p and fld.lSp (note that P{t) (3{t) = by fH.lSp and fl3.2p h 

Notice that 


(4.18) 1 — X < exp(—x), V X e [0,1]. 

(To see this, let W{x) := 1 — x — exp(—x), x G [0,1], then W\x) = —1 + exp(—x) < 0. Hence 
W{x) < 1T(0) = 0 .) 

Applying fl4.18p with x = {j* < j < i), we obtain 

[pj, n] < 1 - p-^ < exp {-pj^) . 

Plugging into fl4.17p yields 

H exp(-p-i)= 2 ’^%xp(- 5 ^ as 

i<j<e \ j*<j<£ / 


with 2 /c + 1 (= K) > P{i), by the well-known result that 


^ oo as £ —)■ cxD . 

i<i<^ 

Since < /ip(£)({/}), we deduce that —)■ 0 as A; ^ cx) , as desired. □ 


4.3. Probability that All Diagonal Entries of the SNF Are 1. 

Theorem 14.71 (along with Theorem 13.6|) implies that the probability that all diagonal entries of 
an SNF are 1 is 0 if m = n; however, as we will see soon, this probability is positive if m < n. 
We will need Theorems 13.21 and 13.81 to determine its value. 


Theorem 4.8. Let E be the n x m diagonal matrix whose diagonal entries are all 1. If m < n, 
then we have 


ri{{E]) = 


n 


i=n—m+l 


C,{^) 




1 , 


if m is hxed 


H 


i=n —m+1 


at 


Y , if 77 , — m is hxed ’ 

7 ' 


as n —)■ cxo 


Proof. Apply Theorem 13.81 with S = {E}, r = m,di = z = l,Sj = 0 for all i,j, and then Theorem 
13.21 with s = 1, ai = m: 


\p,n — m\ 

p P ^ ^ P i=n—m-\-l i=n—m-\-l p 


n n (i-p‘)= n iki-p*) 


n 


i=n—m-\-l 


C(^)’ 


on the strength ofn — m-|-l>2 and the Euler product formula (14.3F 

Finally, thanks to the fact that C(0 i 1 i —)■ cxd , we obtain the limits of p{{E}) as desired. □ 

































SMITH NORMAL FORM DISTRIBUTION OF A RANDOM INTEGER MATRIX 


19 


4.4. Probability that At Most ^ Diagonal Entries of the SNF Are Not 1. 

In this section, we assume that m = n. We provide a formula for the probability that an SNF 
has at most i diagonal entries not equal to 1 and a formula for the limit of this probability as 
n —)■ oo. In particular, when i = 1, this limit is the reciprocal of a product of values of the 
Riemann zeta function at positive integers and equals 0.846936. For bigger i, we prove that this 
limit converges to 1 as £ —)■ cx3 and hnd its asymptotics (see fl4.38p i. 

4.4.1. CycUc SNFs (i=l). 

We shall say that an SNF is cyclic if it has at most one diagonal entry not equal to 1, i.e., if 
the corresponding cokernel is cyclic. Denote the set of n x n cyclic SNFs by 7^ . We will compute 
the probability /i(7^) of having a cyclic SNF, and show that this probability strictly decreases to 
0.846936 ■ • • as n —)■ cxD . As mentioned above, this result was hrst obtained by Ekedahl m Section 
3]. Later Nguyen and Shparlinski [T21 (1-2)] showed that if take a subgroup of Z” uniformly among 
all subgroups of index at most V and let R —)■ cx3, then the probability that the quotient group is 
cyclic is also fi{Tn)- This result is equivalent to computing the probability that an n x n integer 
matrix has a cyclic cokernel using a certain probability distribution different from /r. We do not 
know a simple reason why these two probability distributions yield the same probability of a cyclic 
cokernel. Perhaps there is a universality result which gives the same conclusion for a wide class of 
probability distributions. 


Theorem 4.9. We have 

(i) 


(4.19) 


h'iTn) = 


nr=2C(*) 




(ii) Zn is strictly decreasing in n ; 

(hi) 


1 90 

^2 = ^ ~ 0.923938; 


IV 


C(4) TT^ 

1 I'm ^ — _ 

™ "“c(6)n::4C(^) 


0.846936. 


Proof, (i) Apply Theorem 13.81 with 5 = , r = n — 1, dj = 2 ; = 1, Sj = 0 for alH, j, and then 

Theorem 13.21 with s = l,ai=n,n — 1, respectively: 


(4.20) ^(Tn) = = n ([P'"i+ 


p ^ [p, n] 


[p, l]2[p,n - 1] 


n 


[p, 1 ] 


[p, 1 ] + 


p ^ [p, n] 

[p, l][p,n- 1] 


nr=2C(o 




p ^(1 — p '^) 


1 — p 


-1 


K= 2 a^) 


A n (^+77+Fs+■ ■ •+ 


P^ po 


PT : 


— Zn 


Here in the fourth equality we used the fact that 


(4.21) 


nfei=nn(.-r)=nni.---)-n 


i=2 p 




by virtue of the Euler product formula fl4.3p . 
(ii) We consider the ratio: 


^n+l TT/, .-(n+l)\ 1+P^+P^ + ---+P 


n(i 


p 


') 


-(n+1) 


1 + P“^ + P~^ +-h p~^ 
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thus it suffices to show 
(4.22) 


(1 _ p-(n+i)) . 1+P '+P ' + ---+P < 1 


1 + p~‘^ + p~‘^ H- Vp~ 

for all p. For ease of notation, we denote p~^ by t throughout this paper, then 

tn+l 


LHS of (Km = (1 - r+1) ■ (^1 + 


l + f^ + t^ + '-' + t 
(iii) When n = 2, it follows from dehnition fl4.19p that 


^ < (1 - r+^) (1 + r+^) = l-t2(^+i) < 1. 


Z2 = l[{l-p ') (1+p ') = ^) = 


C(4)- 


(iv) Now assume that n > 3. According to the dehnition fl4.19p of , it suffices to prove that 


n / 1 1 1 

( 1 H- 7^ H- n + ■ ■ ■ H- 

\ p'J p'^ 


C(2)C(3) 

C(6) 


In fact, we will show that 

C(2)C(3) 


(4.23) 


C(6) 


n(i+^+^+ 


hmTT(l + ^ + ^H- 

X X \ pZ pO pTl f 


P P 

We adopt the notation t := p~^. For the left equality of (I4.23p . we observe that 


(4.24) 


l + t^ + t^ + -- - = l + 




1 -t + r 


i + t^ 




l-t l-t ( 1 -^ 2 ) (1-^3) 

Taking the product of this equation over all reciprocals t of primes and applying the Euler product 
formula fl4.3p yields the desire equality. 

For the right equality of fl4.23p . since 
1 + H- \-f- 


° ^ ^ 1 + t 2 + t3 + . . . 

combining with fl4.3p . we obtain 


r+i + r+2 +... r+i + r+2 +... 


i>n 


i + t‘^+ t^+ + 

i + t^ + t^ ^— 


> 


1 + t 2 + t3 - 


n (1 - 


< 


^2 + ^3 . 


= t 


n—1 


^ 1 , as n —)■ cx) 


and complete the proof, where represents a product over all reciprocals t of primes. 

One can also show the right equality of fl4.23p using the fact that 

(4.25) 1 < 1 + p~‘^ + p~^ + ■ • • + p~"' t 1 + P~‘^ + P~^ + • • • , as n ^ cx) 

and the following version of monotone convergence theorem (which will also be very useful later 
in proving Theorem I4.13l fiiii). 

Theorem 4.10. If real numbers Xij {i,j = 1, 2,...) satisfy 1 < Xtj f Xi as j ^ oo for all i, then 
we have 


(4.26) 




7—>-CxD 

i=l 


i=l 


Here we allow the products and the limit to be infinity. 
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Proof. Applying the monotone convergence theorem to logXjj (> 0) gives 


CXD OO 



Thus 

OO OO 

lim log IT Xij = log IT Xi , 
7^00 

i=\ i=l 

and fl4.26p follows. 


□ 


Thanks to fl4.25l) . we can apply Theorem 14.101 with Xjj = + and Xi = 

1 + + • • •, and arrive at the desire equality. □ 

Remark 4.11. (1) The proof of Theorem 14. 9l fiv) is reminiscent of (though not directly related to) 
pT| Exercise 1.186(c)]. 

(2) Theorem I4.9l (i). (iv) and the numerical value of (hi) are obtained in [HI Section 3] via a 
slightly different approach. We have provided a complete and more detailed proof. 


4.4.2. More Generators (General i). 

Now we consider the SNFs with at most i{< n) diagonal entries not equal to 1, i.e., whose 
corresponding cokernel has at most i generators. Denote the set of such n x n SNFs by TniG). In 
particular, when i = n, we have niTniji)) = 1. The above discussion on cyclic SNFs is for the 
case £ = 1. We will compute the density p(7^(£)) and its limit as n ^ cx), show that this limit 
increases to 1 as £ —)■ cx), and establish its asymptotics. 

We start with a lemma which will play an important role in our proof (as well as in Section 15.21 
below). 

Lemma 4.12. For any positive number x < 1/2, the positive sequence {[1 /x, is decreasing 

and thus has a limit as k ^ oo: 


(4.27) C{x) := (1 - t) (1 - x^) • • • G l) . 

This also implies that C{x) ^ 1 as x ^ 0 and that [1/x, k] G 1) for all x E (0,1/2] and 

k>l. 

In particular, when x = 1/p, we have 

(4.28) [p,k] fCp ■.= C{l/p) E \e~‘^^^^~^\l) C |]e“^, l) , as k ^ oo , 

Cp —)■ 1 as p ^ OO , and \p, k] G 1) for all p and k > 1. 

Proof. The sequence [1/x, k] is strictly decreasing in k because 0 < 1 — < 1 for all j > 1. 

To get the lower bound for C (x), we will use the following inequality: 

1 -y 


(4.29) 


hip > 


y 


, VpG(0,l]. 


(To see this, let tp{y) := Inp + (1 — y)/y, then 'tp'{y) = 1/y — 1/y^ < 0. Hence i/(p) < -0(1) = 0.) 
Applying (I4.29p with p = 1 — x-^ (j > 1) yields 


(4.30) 


In (l — x-^) > —Y 


x-' 


x^ 


> - 2 x^ 


as x-^ < 1/2. Summing up fl4.30l) over j from 1 to A:, we get 

k OO 

In [1/x, k]> — ^ 2 x^ > — ^ 2 x^ 
i=i i=i 

Hence C(x) = \imk^oo[l/x, k] > 


2 x 


1 — X 


□ 
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Theorem 4.13. We have 

ii) 


(4.31) 

where 


(4.32) 


h^{Tnm = n Zn{p. 1) = n ’ 

„ lh=2SGj „ 


Zn{p, i) = Pp{Tn{i)) = \p,n]Y^ 


p [p, n] 


[p,i]‘^[p,n-i\ ’ 


i=0 


Yn(p, t) = Whz^p, () = [p, 1] 5^ WW 

\pw] 


p * [p, n] 
n — i] 


[P, 1 ] 


(4.33) F„(p,£) t [P, 1 ]X^ 7 — r^=-y{.pJ) as n^cx), F(p,£)t^ 

i=o 

and 

(4.34) 


as L ^ oo , 


a 


£ _j2 

p 


Pp{Tn{()) = Zn{p,tj^ Y-^Y{p,tj=Cp^- — —=:Z{p,tj as n ^ oo , 

\p,i] ~^\P^A 

where Cp = (1 — p“^)(l — p~^) ■■■ as defined in fl4.28p and fl4.27p . then it follows from fl4.33p that 
(4.35) Z{p,i)'\'l as £ —)■ oo ; 

(iii) 

1 


(4.36) 


piTnii)) = Zn{i) ^ 


n y(p, 0=n ^(p- <■) =■ 2(«) 

lh=2SGj „ „ 


as n —)■ oo , 


and Z{t) f 1 as i ^ oo ; 
(iv) 


(4.37) lim PpiXii)) = Z{p, i) = l- 

n—>-cxD ^ 

more precisely, this O (p~^^) G (0, 2 p~‘^^) ; 


p^ — p 


, p-^ + O (p- 2 ^) 


as i ^ oo ; 


V 


(4.38) 


lim p(Tn{i)) = Z(i) = 1 - n -2-^ + 0 (4-^)1 as £ ^ oo , 

' L \ / J 


where ^ 3.46275 . 

Parts (ii) and (iv) also hold with p = 1/x for any x G (0,1/2]. 

Figure[T]and Table[T]below illustrate the asymptotics (14.381) of Z{£) and fast rate of convergence. 

Remark 4.14. The convergence result fl4.35p in Theorem I4.13l fii) with p = 1/x implies Euler’s 
identity: 

(1 — xfi{l — x'^fi • • • (1 — x^fi (1 — a;)(l — x'^) ■ ■ ■ 
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Figure 1. Asymptotics of Z{tj 


> 

3.6 

;2h+W(l - Z(£)) 

-ln(l -C 2 (1-Z(£))) /ln2 
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Table 1. Asymptotics of Z{t) 


i 

Z{i) 

1 - Z{i) 

Z{i)) 

-ln[l - C '2 2 (^+^)'( 1 -Z(£))] /ln 2 

1 

0.846935901735 

1.53064098265 x 10-^ 

2.44902557224 

1.77225611430 

2 

0.994626883543 

5.37311645734 x lO-^ 

2.75103562616 

2.28255339912 

3 

0.999953295075 

4.67049248389 x 10"^ 

3.06085395424 

3.10703467197 

4 

0.999999903035 

9.69645493161 x 10"® 

3.25359037644 

4.04926385851 

5 

0.999999999951 

4.88413458245 x 10-^^ 

3.35635172814 

5.02441603986 

6 

1.000000000000 

6.05577286766 x lO"^® 

3.40909705378 

6.01220652280 

7 

1.000000000000 

1.86255532064 x lO-^^ 

3.43580813230 

7.00610418193 

8 

1.000000000000 

1.42657588960 x lO-^^ 

3.44924885316 

8.00305233425 

9 

1.000000000000 

2.72629586798 x lO-^o 

3.45599059345 

9.00152622794 

10 

1.000000000000 

1.30126916909 x lO’^® 

3.45936681921 

10.0007631292 


Proof, (i) The first equality follows from Theorem 13.81 with S = Tn{i), r = n — i, di = z = 1, 
Sj = 0 for all i, j, and Theorem 13.21 with s = 1, oi = n,n — 1,... ,n — i, respectively. 

The second equality follows from dehnition fl4.32l) and fl4.21l) . 

(ii) We observe that 

, , = (1 - p-^) (1 - ■ ■ ■ (1 - I 1 as n^oo. 

This leads to the hrst result of fl4.33p . 

Since W(p, P) is also increasing in I by dehnition (I4.32p . so is Y{p, i), and for all i < n, we have 
(4.39) Ye{p,i)<Yn{p,i)<Y4p,n). 

Further, we derive from 

1 = PpiTnin)) = ^^W(p,n), 

that 

Yn(p,n) = 1^’ and similarly, YAp,i) = . 
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Plugging into fl4.39p . we obtain 


Taking n —)• cx) yields 


Ip-1) < < M ^ [p. 1] 




[P, 1 ] 


<Y{pJ)< 


[p, n] Cp 

[P, 1 ] 


a 


[p, i] 

Then taking £ —)■ cx) and applying Lemma 14.121 leads to the second result of fl4.33p . 

Finally, on the strength of fl4.33p and Lemma [4. 121 we obtain fl4.34p from dehnition fl4.32p : 

Zn{p,i) = j^Yn{p,i) ^ j^Y{p,i) as n^oo. 

IP, 1] [P, 1] 

This proof also carries over to p = 1/x for any x G (0,1/2], 

(iii) It follows from dehnitions fl4.3ip and fl4.32p that 

[p,n 


(4.40) 


znV) =n ^»(p 4)=If phi n ^“(P’ 


By virtue of fl4.2ip . we have 

|p.n] 


(4.41) 


n 


1 1 


[p,i] nr=2C(o E2C(*) 


0.435757 as n —)■ cxo . 


Further, this limit 


_i OO OO 

TEcw ^ i n (1 - P-)=n n (1 - P-)=n yij 


Hence 

(4.42) 


n [Pi TT 

^ [P, 1] [P, 1] 


as n —)■ OO . 


since 


We can also deduce fl4.42p from Theorem 14. 101 with xtj = [pi, l]/\pi,j] si 

1 / [Pi 1] ^ [Pi 1] 

1 < 7-7 j —— as n ^ OO 

[p, nj Cp 

by Lemma [4. 121 

For the second product on the right-hand side of (14.401) . from fl4.20p in the proof of Theorem 
I4.9l (i). we see that W(p, 1) = 1 +p“^ +p~^ + ■ ■ ■ +p“” > 1. Since W(p, i) is increasing in i, we have 
Yn{p,i) > 1 as well. In conjunction with (I4.33p . we can apply Theorem 14.101 with Xij = Yj{pi,i) 
to obtain 

(4.43) WYn{p,tj "[WY{p,tj as n ^ oo . 

p p 

Plugging fl4.4ip . (14.431) and (I4.42p into (I4.40p along with dehnition (I4.34p yields (I4.36p : 

(4.44) z„(<) ^ Y(p, <)=n 7 % n ^(p- ^)=n 7 %i ^(p- <■)=n ^(p- 

n *=2 c(*) V [p, 1 ] V V [p, 1 ] V 


as n —>■ cx). 
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Since Yn{p,i) > 1 and Yn{p,i) is increasing in £, so is Y{p,i) (recall fl4.33p l. Thus we can apply 
Theorem 14.101 with Xij = Y{pi,j) to obtain 

Yly{pj) t n ^^ oo. 

Op 

p p 

Finally, we plug this into the second expression of the limit of Zn{i) in fl4.44p : 

z{i)= iimz„(£) = rr^rry(p,£) 11 as i^oo. 

n—>-oo W, 1 

p 1^ ‘ -I p 


(iv) We prove for the more general case p = 1/x with x G (0,1/2], Let 

Recall that C{x) = (1 — x)(l — x^) • • • and [1/x, i] = (1 — x)(l — x^) • • • (1 — x*). 
Since R(x, £) = Z(l/x, £) 11 as £ —>■ oo by fl4.35p . we have 


OO ,-2 

X 


^ M P2 


X 


OO ,-2 

X 


X _ V (x, i) 

\'\ I nr i\2 C^A Y\ 


OO ,-2 

X 


C{x) ^[1/T,i]2 ^[l/a:,i ]2 C{x) ,^Jl/x,i] 


Thus for any x G (0,1/2], we obtain 

x-(^^i)^C'(x) [1 — R(x, £)] = x“^^'''^^^C^(x) 


1 V{x,i) 
“ C{x) 


- rr-A+^P C‘^ix)x^ 
^ [l/x,i ]2 


= X 


i=i+\ 


OO OO 

>2 


^ x^ n (1 “ = n 


OO OO 

2 „•2_(•/l_l_1^2 


X 


F-p+ 1 )^ 


n 


*=£+! j='i+l j=£+2 i=£+2 j=i+l 

- X 2x^+2 


1 — 2 ^ ^ x^ + Ai I + A 2 — 1 — 

i=£+2 


1 — X 


+ A 1 +A 2 , 


(4.45) = 
where 

(4.46) 0 < A 2 := JJ (l - x^')" < ^ 

and 


OO OO 

,-2_/'* I 1 'i2 


2£+3 


i=£+2 


j=i+l 


OO OO 

2 , V- ^P-A+i? ^ J2 Y = - - < x^^ 

i=i+2 i=2e+3 


1 — X 


(4.47) 0 < Ai := JJ (l - x^Y - ( 1 - 2 ^ xM < 4 ^ 


xi+f = 


j=e+2 


j=i+2 


j,j'>e +2 


4x2r+4 

(l-x )2 




as 0 <x<l/ 2 , thanks to the inequality: 


0<n(i-U- (i-E'*') s E ■’A 


2=1 


2=1 


for (5i, (^ 2 ,..., G [0,1], which can be proved easily by induction on u (the left inequality was 
proved in fl2.18p . For the right inequality, base cases: m = 1,2; inductive step from m to m + 1: 

(1 - (5„+i) nr=i(i “ ^ (1 “ <^«+i)(i - DLi 1 - 

^u+l A 1 — X]i=l + 'I2l<i<j<u+l^i^j)- 

Plugging fl4.46p and fl4.47l) into fl4.45l) yields fl4.37l) . 





































26 


YINGHUI WANG AND RIGHARD R STANLEY 


(v) Since Z(£) = flp Z(p, i) by definition fl4.36p and 0 < Z{p,i) < 1 for all p , we have Z{£) < 
Z{2,£). Thus it follows from (iv) that 

(4.48) Z{£) < Z(2, £) = l- [l -2-^ + 0 (4“^)] as £ ^ oo . 

On the other hand, we notice that when £ >2, the 0{p~‘^^) in (14.371) satisfies 

0{p-^^) < 2p-^^ < 4 • • P ~^) 

pZ p^ — p 

thus 

Z(p,i)>l-Cpp-f‘*'>\ 

Hence 

(4.49) Z(e) = n Z(p. 0 > Z(2, <) n (i - G“V“<'+‘’’) > 1 - (1 - Z(2, Q) - Cpp-''‘+»\ 

P P>3 p>3 


Here we took advantage of the inequality fl2.18l) . Thanks to fl4.28p . the positive sum 


p>3 


-m-p+i)2 


p>3 


p>3 




p 


{i+i? 


p>3 


2 \ ^ V 2 ^ ^ 


p 


( \ t 

■4 = 2-('+')^O(4-0 

Finally, combining with fl4.48p and fl4.49p leads to fl4.38p . 


□ 


Remark 4.15. When £ = 1, in the proof of Theorem 14.91 we wrote y(l/x, 1) as (1 — a:®)/(l — a:^)(l — 
x^) (see (14.2411 ) in order to represent Z(l) = ^Xp^iPt 1) /YW= 2 C{P) as the reciprocal of a product 
of values of the Riemann zeta function at positive integers. However, this is not the case when 
£ > 1; in fact, in general Y{l/x,£) is not even a symmetric function in x, for instance. 


Y 



1 — X — + 2x^ — X® + X® 

(1 — x)^(l + x )2 


Y 



1 — X — x^ + 2x^ + X® — 2x® — x’^ + X® + X® — x^^ + x^^ 
(1 — x)®(l + x)2(l + X + x'^y 


5. Properties of the SNF Distribution Function pps 

In this section, we first fix p, s, m, n and find the maximum and minimum of the probability 
density function pps of (13.11) . Then we free p,s,m,n and study the monotonicity properties and 
limiting behaviors of Pps[{Da}), as a function of p,s,m,n and a (recall from Theorem 13.21 the 
notation of vector a = (oi, 02 ,..., a^) as well as its corresponding diagonal matrix Da € §). 

For convenience, we replace m — a, by fej (0 < i < s) in (13.111 to get a simpler expression for 

Rp^iiDa})- 


(5.1) 


f{p,s,m,n',b) :=p 


Yi=iD'+bi)bi 


\p, n' + m] \p, m] 

[p, n' + bs] [p, bs] riLi [P) bi-i - bi] ■ 


Here and throughout this section, we shall assume that p is a prime, that s, m and n are positive 
integers, that n > n' := n — m > 0, and that b := ( 61 , 62 , • • •, bg) is an integer vector satisfying 
TTT' = bo > bi > ■ ■ ■ > bg > 0. 
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5.1. The Maximum and Minimum. 

We show that f{p, s, m, n', •) attains its maximum at either (0,0,..., 0) or (1,1,..., 1) depending 
on p, s, m and n', and its minimum at (m, m,, m). 

Theorem 5.1. For fixed p,m,n and s, the maximum and minimum of f{p, s,m,n', ■) are given 
as follows. 

(i) If p > 2 , s > 1 or n' > 0 , then 

, , . [p,n' + m] 

max /(p, s, m, n , b) = —^--— , 

b [p, n'J 

and the maximum is achieved if and only if b = ( 0 , 0 ,..., 0 ) := 0 , in other words, if the corre¬ 
sponding matrix Da is full rank; 

(ii) If p = 2 , s = l, n' = 0 and m > 1, then 

fi ' u\ [2,mf 

mr ^ [2.1] [2. „-!] ■ 

and the maximum is achieved if and only if b = ( 1 ); 

(iii) In both Case (i) and Case (ii), we have 

min f{p, s, m, n', b) = p”**-"' ^ 

b 

and the minimum is achieved if and only if b = (m, m,..., m), in other words, if the corresponding 
matrix Da is the zero matrix. 

(iv) Ifp = 2,s = l,n' = 0 and m = 1, then 6 = ( 1 ) or ( 0 ), and they have the same value of 
/■• 1 / 2 . 

Proof, (i) We proceed by the following two lemmas which show that the bfs are all equal at the 
maximum of /(p, s, m, n', •), and that 6 * = 0 or 1 depending on p, s, m and n'. 

Let b = (5i, 62 , • • •, bg) be an arbitrary s-tuple with m = h^ >hi> ■ ■ ■ >hs>^ . 


Lemma 5.2. If bi > bj+i for some i G {1, 2,..., s — 1} (s > 2), then we have 

f{p,s,m,n',b') > f{p,s,m,n',b), 

where b' = {b[,b 2 ,... ,b'g) with 6' = 6, — 1 and 6) = bj for all j i. Note that b' still satisfies 
m = Fq > b[ >■■■> b'^ > 0 . 

Lemma 5.3. Let ip{b) := f{j), s, m, n', {b, b,..., b)), 0 < b < m, then for all 0 < b < m, we have 

<1, ifp = 2,s = l,n' = 0,m>l and 6 = 0 
= 1, ifp = 2, s = l, n' = 0, m = l and 6 = 0. 

> 1, otherwise 


Pjb) 

(p(6+ 1) 


These lemmas are proved right below this proof. Thanks to Lemma 15.21 the maximum point of 
/(p, s, m, n', •) must have the form (6, 6,..., 6) with 0 < 6 < m. Therefore it reduces to hnding 
the maximum of (p(-). 

Since p > 2, s > 1 or n' > 0 , it follows from Lemma 15.31 that 


(5.2) (p(0) > <p(l) > • • • > <p(m). 

Hence the maximum of (p(-) is (p(0) = , as desired. 

(ii) When p = 2,s = l,n' = 0 and m > 1, it follows from Lemma 15.31 that 

(5.3) (p(0) < (p(l) and (p(l) > • • • > p{m). 
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Hence the maximum of (p(-) is (^(1) = [2 desired. 

(iii) We proceed by the following lemma (proved right below this proof) which shows that at 
the minimum of f{p, s,m,n', •), all the fcj’s {i > 1) equal m. 


Lemma 5.4. If bi < hi-i for some i G {1, 2,..., s — 1} (s > 2), then we have 

f{p,s,m,n\b') < f{p,s,m,n\b), 

where b' = {b[,b 2 ,... ,b'f) with 6' = 6* + 1 and 6' = bj for all j ^ i- Note that b' still satisfies 
m = 6g > > ■ • • > > 0 . 

Thanks to Lemma 15^ the minimum point of f{p, s, m, n', •) must have the form [m, m,... ,m,b) 
with 0 < 6 < m . Further, since f{p, s, m, n', (m, m,... ,m,b)'j = f|5.ip L 

where p is dehned in Lemma [5.31 with s = 1, it reduces to hnding the minimum of ^{■)- 

Case (i) When p > 2 or n' > 0, it follows from fl5.2p that the minimum of (p(-) is (p(m) = 

p—{n'-\-m)m 

Case (ii) When p = 2, n' = 0 and m > 1, it follows from fl5.3p that the minimum of (p(-) is 
min {(p(0), (p(m)}. Since 

(p(0) = [2, m] > (1 - 2"^)”' = 2-^ > 2"™' = p{m) , 
the minimum of (p(-) is still ip{m). 

Hence the minimum of / is always p-^C+m)m achieved at (m, m,... ,m). □ 


Proof of Lemma 15.M It follows from dehnition flS.ip that 

f{p,s,m,n',b') ^ _ [p, bi_i - bj] [p, bj - k+i] 


> p ■ 


f{p,s,m,n',b) ^ 

[p, bi i - bi] [p, bi - bi+i] 


[p, bi_i - 6'] [p, 6' - 6i+i] 


[p, bi_i -bi + 1] [p, bi-1- bi+i] 


= P - 


p 






>p(l—p ^)=p — 1>1, 


as desired, where in the second last inequality, we used the condition that bi > 6j+i to get 1 — 

-(bi-bi+i) y I _ p-1 _ □ 


p 


Proof of Lemma \5.3[ By the dehnition of <p and fIS.ip . we obtain 

pjb) ^ ^s[{n'+b+i){b+i)-{n'+b)b] . [p, n' + b + l][p,b + l][p,m - b - 1] 

(p{b + 1) [p, n' + 6] [p, 6] [p, m — 6] 

(5.4) = pd^'+2b+i) . {^-P ^^ ^ ^ > p^C'+2b+i) _ p- 1^2 ^ 

where we used the fact that 

1 _ p-{n'+b+i)^ ^ _ p-ib+i) > 1 _ p-i and 1 - pCm-b) < 1 _ 

Case 1. s{n' + 2b + 1) > 2 . 

The right-hand side of fl5.4p is at least 

p^ (1 — p“^) ^ = (p — 1 )^ > 1 . 

Case 2. p > 3 . 

The right-hand side of (15.4p is at least 

p (1 - p-^Y > 3 (1 - 3"^)^ = 4/3 > 1. 

Case 3. p = 2 and s{n' -|- 26 -|- 1) = 1, which requires that s = 1 and n' = 6 = 0. 
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Plugging into fl5.4p yields 


(^( 6 ) _2{l-2-^f _ 1 /<!, ifm>l 

^{b + l)~ 1 - 2 -™ ifm = l 

and completes the proof. 


□ 


Proof of Lemma By definition fIS.ip . we obtain 

f{p,s,m,n',b') ^ _ [p, bi_i - bj] [p, bj - bj+i] 

f{p, s, m, n\ b) [p, bi_i - 6'] [p, b\ - h+i] 

< „-i . [p,bi-i-bi][p, bi-bi+^] ^ _ 1 _p-(fei-i-bd ^ ^ 1 _ ^ ^ 

“ [p, bi-i -bi-1] [p, h + l- bi+i] 1 - p-ih+i-h+i) 1 - p-i p - 1 “ 

as desired, where in the second last inequality, we used the condition that > 6j+i to get 1 — 

p-(bi+l-bi+i) y I _ p-1 _ □ 

5.2. Monotonicity Properties and Limiting Behaviors. 

Now we free p, s, m and n'. We will see that the monotonicity properties and limiting behaviors 
of / of fl5.1l) when b = 0 (i.e., the corresponding matrix Da is full rank) differ tremendously 
from those when b ^ 0. Specihcally, we show that / is increasing in n',p and decreasing in m 
when 6 = 0 iTheorem 15.5p . but decreasing in n' and increasing in m when b ^ 0 fTheorem 15.6p . 
Further, with regard to limiting behaviors, when 6 = 0, the limit of / as p, m or n' —)■ cxo is 
positive (note that / is independent of s) fTheorem 15.51) : whereas when 6 7 ^ 0, the limit of / 
is still positive as m —)■ 00 or s —)■ cx 3 with X]i=i bi bounded (Theorems 15.10115.11|) . but zero as 
max {p, n', h} ^ 00 (Theorems 15.7115.91) . Lemma 14.121 is crucial in the analysis of limiting 
behaviors of /. 


5.2.1. The Case of b = 0 . 
Let 


We derive the following monotonicity properties and limiting behaviors of /o with the help of 
Lemma 14.121 


(1 -p ^) . 


(5.5) 


/o(p,m,n') := /(p, s, m, n', 0) = 


[p, n' + m] 
[p,n'] 


n'+m 

n 

j=n'+ 


Theorem 5.5. The function fo{p,m,n') of fl5.5p is strictly increasing in p,n' while strictly de¬ 
creasing in m, and satisfies 

C 

lim fQ{p,m,n') = -—^ < 1 , lim inf/o(p, m, n') = 1 and lim inf fQ{p,m,n') = 1 , 

m—>-00 P, TlO n'—>00 m p—>oom,n' 


where Cp = {1 — p ^)(1—p ■ as defined in Lemma [4-1^ In particular, we have 


lim /o(p, m, 0 ) = (Fp and 


lim /o(p,m,n') = 1 = lim /o(p,m,n') ; 
n'—>-oo p—>-00 


the first eguality characterizes Cp as the limit of the probability that a random mxm integer matrix 
over 'Ljp^'L is nonsingular as m ^ 00 . 


Proof. Utilizing the expression on the right-hand side of fl5.5p . we obtain the monotonicities. Then 
we apply Lemma 14.121 to get 

C C 

inf /o(p, m, n') = lim /o(p, m, n') = — ^ ^ = 1 as n' ^ cx), 

m m^oo n J Cp 
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and 


inf fo{p, m, n') = lim /o(p, m, 0 ) = Cp —)■ 1 as p —)• cx). 

m,n' m—>-oo 


□ 


5.2.2. The Case of b ^ 0 . 

We first present the monotonicity properties of f{p, s, m, n', b) in n' and m. 

Theorem 5.6. Suppose that 6 7 ^ 0 . The function f{p, s, m, n', b) is strictly decreasing in n' while 
strictly increasing in m. 

Proof. Recall that hi>h 2 >---hs>t]. Since h 7 ^ 0, we have 61 > 1. Thns the ratio 

/(p,g,m,n^ + l,b) ^ _ \p,n' + l + m][p,n' + bf\ 


= P 


f{p,s,m,nfb) ^ 

_p—(n'+l+m) 


[p, n' + m] [p, n' + 1 + 5^] 


= n . 


< p 


-1 


< 1 , 


and 


1 _ p-{n'+l+bs) ' ^ I _ p-1 p — I 

f{p, s,m+ 1, n', b) \p,n' + m + l][p, m + l][p, m — bi] 

f{p, s, m, nf b) [p, n' + m] [p, m] [p, m + 1 — bi] 

_ p-in'+l+m)^ _ p-(m+l)^ _ p-(m+l)j 2 ^ ^ ^p-im+l) 


> 


1 — p—(m+l—bi) X — p“"^ 1 — p“ 

as p > 2 . 

Recall from dehnition fl5.ip that / is the prodnct of a power of p 

fi{p, s,m,n',b) : = p-T,i=iC+bi)bi 

and a fraction 

[p, n' + m] [p, m] 


> 1 


□ 


(5.6) 


f 2 {p,s,m,n',b) := 


[p, n' + bs] [p, b,] n ■=! [p, bi-i - bi] ' 

When s is hxed, thanks to Lemma 14.121 the fnnction /2 dehned in 05.61) is bonnded regardless 
of the valnes of other variables. Moreover, when m (instead of s) is hxed, this resnlt also holds 
since — h) = m — bg < m implies that 

]^[p,6i_i - 6i] = ]^ JJ ( 1 -P"^) > ( 1 > 2"”^. 

2=1 2=1 j = l 

These observations lead to the following zero limiting probabilities. 

Theorem 5.7. We have 

max f{p, s, m, n', b) = 0 when s is fixed 


and 


lim 

max {P ,n', 


lim /(p, s, m, n\ 6 ) = 0 when m is fixed. 

max {p. n', 6i}-)-oo, 670 


Proof. When s or m is hxed, we have shown that /2 is bonnded. On the other hand, we have 
fi{p, s,m,n',b) = p~'Pi=iC+bi)bi _ p-«'Ei=Oi-Ei=Oi 
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as long as b 7^ 0 and 

( 5 . 7 ) 

Noticing that 


max < p, n' bf 


—)• oo. 


i=l 2=1 


^ bi < n' ^ bi + ^ 6^ < n' ^ bi + [ y^hj 


2=1 


2 = 1 2 = 1 

/ 


2=1 


. 2 = 1 


thus fl 5 . 7 p is equivalent to max {p, n', X]i=i bi} —>■ oo . 


□ 


Remark 5 . 8 . Let r (< s) be the number of nonzeroes in {6i, 62, • • •, i.e., 6^ > 0 = hr+\ (we 

dehne 6^+1 = 0 ), then rhr,bi < X]i=i decreasing property of the bj’s. Hence 

bi —t 00 if and only if max {bi, r} —>■ 00 . In particular, when s is fixed, we have X]i=i b* —)■ 00 
if and only if 61 ^00. 


Moreover, if we free s, m, n' but fix p and let b* —)• 00 , then / also goes to 0 . 


Theorem 5.9. For a fixed prime p, we have 

lim max f{p, s, m, n', 6) = 0 . 

ELi m,n' 

Proof. Since X]i=i b* —)■ 00 , we can assume that b 7 ^ 0. Moreover, if r (< s) is the number of 
nonzeroes in {bi,b 25 • • • then max { 61 , r} —)■ 00 (see Remark [5.8p . which is equivalent to that 
bi —)■ 00 or r ^ 00 holds. 

Case 1. bi —>■ 00 . 

For any fixed p, s, m and n, from Lemma [5.21 we see that for h' = (bi, bs, bg,..., bg), 

/(p, s, m, n', b) < /(p, s, m, n', b') 

= p-(n'+ 6 i)bi-(s-i)(n'+b.)b._ [p,n^ + m][p,m] _^ _ 1 

[p, n' + bg] [p, bg] [p, m - bfi [p, bi - b*] “ 
on the strength of Lemma 14.121 Hence 

X 

max /(p, s, m, n', b) < p“^^ • -—^ t 0 , as bi ^ cxd . 

m,n' (e 

Case 2. r ^ 00 . 

For any hxed p, s, m and n, from Lemma 15.21 we see that for h' = {br, br,..., br, bg, bg,..., bg) 
(with r bfis and (s — r) b^’s), 

/(p, s, m, n', b) < f{p, s, m, n', b') 

= r^-r(n'+br)hr-(s-r)(n'+hs)bs _ [p, u' + m] [p, m] _ ^ _ 1 

[p, n' + b^] [p, b^] [p, m - br] [p, br - b^] “ ’ 

on the strength of Lemma 14.121 Hence 

max fip, s, m, n', b) < p~'' ■ - — — —?• 0 , as r —)■ oo . 

m,n' (e 

□ 


All the limits of / we have found so far equal zero. To attain a nonzero limit, we must have a 
bounded max {p, n', h}- We may £x p, s, n', b, let m ^ oo and apply Lemma 14 . 121 . 

Theorem 5.10. For fixed p, s,n' and b 7 ^ 0 , we have 


lim f{p, s, m, n', b) = p 'Fi=iC+bi)bi 

m^oo 


_ Cp _ 

[p, n' + b,] [p, b,] 11^2 bi_i - bi] ‘ 
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We may also weaken the constraints by fixing p, n' and only. A natural way to achieve 

this is to fix the first few say bi,b 2 ,... ,br (r < s fixed), and set the rest to be zero no matter 
how big s is. According the definition ( 15 .Ih of /, for b = ( 6 i, 62 , • ■ ■, 0, 0,..., 0), we have 

(5.8) /(p. ra, t) = p" , . 

[P^n'] Oi-1 - bi\ 

which is independent of s. Coupling with Theorem 15.91 gives the following. 


Theorem 5.11. When m,n' and p are fixed, for any given infinite integer sequence {bo,bi, 
with m = bo > bi >■■■> bi > 6^+1 > • • • > 0 , we have 


lim f{p,s,m,n',b^) = 


0 , 


p 


-U=W+bi)bi . 


[ p , n '+ m \ [ p , m \ 


if bi^oo 
otherwise 


} 


where b^ := {bo, 61 ,..., bs) and in the second case, r is the number of nonzeroes in { 60 , &i, • • • } and 
finite (see Remark lST^) . and A+i = 0 . 
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